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Abstract 

The large-n expansion is applied to the calculation of thermal critical exponents describing the 
critical behavior of spatially anisotropic d-dimensional systems at m-axial Lifshitz points. We 
derive the leading non-trivial 1/n correction for the perpendicular correlation- length exponent 
VL2 and hence several related thermal exponents to order 0(l/n). The results are consistent 
with known large-n expansions for d-dimensional critical points and isotropic Lifshitz points, as 
well as with the second-order epsilon expansion about the upper critical dimension d* = 4 + m/2 
for generic m £ [0, d]. Analytical results are given for the special case d = 4, m = 1. For uniaxial 
Lifshitz points in three dimensions, 1/n coefficients are calculated numerically. The estimates of 
critical exponents at d = 3, m = 1 and n = 3 are discussed. 
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1. Introduction 



The (multi)critical behavior of strongly anisotropic systems in the vicinity of a Lifshitz 
point (LP) [1,2] is given by the Landau-Ginzburg- Wilson (LGW) Hamiltonian 



H[<t>] = J 



d d .i 



4! 



(1.1) 



where = 4>(x) is a classical n- vector order-parameter field, and the interaction is given 
by a standard 0{n) symmetric 4 term. 
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Important is that physical properties along different spatial directions essentially differ 
and this difference cannot be removed by simple rescalings of the theory. 1 Accordingly, 
the d-dimensional coordinate space is split into two Euclidean subspaces R d ~ m and K m . 
Each position vector x e R d ~ m x R m has a (d — m)-dimensional "perpendicular" com- 
ponent x± and an m-dimcnsional "parallel" one, x\\. The gradient operators V_l and 
V|| act in M. d ~ m and M m , respectively, and A|| = V| is the Laplacian associated with 
parallel directions. The term (A^cf)) 2 is fully rotationally invariant in K m . Generically, 
less symmetric fourth-order derivative terms are allowed to contribute into (1.1) thus 
introducing an additional anisotropy in M. m if m > 1 [8] . These can be written as a linear 
combination w a r i " fc j(V|| ;i V|| j: ,</))(V|| ) fcV|| j /(^) where T are totally symmetric fourth-rank 
tensors compatible with the symmetry of the system. 

When the number of anisotropy axes m shrinks to zero we retrieve an isotropic § d 
theory with the upper critical dimension d* = 4 and the usual critical-point (CP) behav- 
ior. Another marginal situation is when m extends to d. The resulting <fim=d theory with 
d* = 8 corresponds to the isotropic LP [1,9,10]. For generic m £ [0,d], there is a line of 
upper critical dimensions d* = d*(m) = 4 + m/2. 

The cf) A coupling uo, as well as the anisotropy parameter do, are positive constants, 
whereas To and po are allowed to change sign. As usual, To is a linear function of temper- 
ature, while po is related to an additional " non-ordering" field specific for the underlying 
physical system (for a review see [11]) and controls the crossover from LP. 

In order to reach the LP, the parameters To and po have to be tuned to their special 
values Tq P and p LP . In the Landau approximation, t lp = p LP = 0. Similarly, in the 
renormalization-group theory the both renormalized counterparts of To and po vanish at 
the LP: t = p = 0. 

It is demanded that the terms (Vj_0) 2 and cr ( A|| 0) 2 scale in the same way. Hence, in 
a cutoff-regularized theory, integrations over perpendicular and parallel momenta k± and 
fey have to be restricted asymmetrically via |fcj_| < A and |fc||| < er^^A 1 / 2 . For large 
A, the shifts r LP and Pq P behave as oc A 2 and cx A, respectively, while in dimensional 
regularization (see e. g. [12,13]) they vanish. 

A simple dimensional analysis (cf. [14-16]) implies that 

[to] = -2M, [x\\] = \[ Xl ] + \[ao], [4>] = -\[d d x] + [xj = I(rf-|-2)[ M ]-|[a ], 

I 771 

[po] = 2[x\\/x ± ] = -[x ± ] + -[<7 ], [u ] = [p e ] + -j[(t ], where [p] = -[x±], 

p, being an arbitrary momentum scale. The deviation from the upper critical dimension 
e is given by e = d*(m) — d = 4 + m/2 — d. The (dimensionful) parameter a^ 4 provides 
the scale for "measuring" distances in the parallel subspace. Its classical dimension is 
given by [cr^ 4 ] = [xn/x 1 / 2 ]. The dimensionless combination v = a^^x^ /x 1 / 2 shows up 
as the argument in the scaling function of the free propagator written in the scaling form 
G^(x) = x^ 2+E a V) m ^ 4 $(u). 2 The effective expansion parameter uo = uou m ^ has the 
usual ^-dimension [p 6 ] and becomes marginal at d = d*(m). Similarly, the scaled field 

1 For discussions on this issue see [3, Ch. 4.2], [4], [5, pp. 15-17], [6, pp. 13-14], and [7] in the context of 
Lorcntz violating theories. 

2 For general m and d, the function has been derived in [15,17]. Its explicit expressions in different 
special cases can be found in [18,19,17,20,21]. 
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4>{x)a™' becomes dimensionless at the lower critical dimension [22,11] dg(m) = 2 + m/2 
appropriate for the models with continuous 0(n) symmetry in the case n > 1. 

Provided that dt{m) < d < d*(m), < m < d, and the LP exists (see [11,20]), the full 
two-point correlation function of the interacting theory obeys anisotropic scaling via 



is the scaling dimension of the field, and all non-universal metric factors are omitted. 
The anomalous field dimension r)L2 is similar to the usual Fisher exponent rj of the pair 
correlation function at CP. The key difference of (1.2) from conventional CP scaling 
forms is in that the distances along perpendicular and parallel directions are rescaled 
with distinct scale factors: Here the anisotropy index 9 is different from 1. We stress also 
that 9 deviates from its classical value 9q = 1/2. The difference, 9\ = 9 — 1/2, is of order 
0(e 2 ) in the epsilon expansion [15] and 0(l/n) in the large-n expansion [20]. Its non-zero 
value is provided by the non-trivial renormalization of the parameter <7 . 

We would like to recall that according to the definition of [1] (see also reviews in 
[2, 23, 11], [24, Ch. 6]) there is the following physical picture of the LP behind the above 
formal description. In underlying systems, two different low-temperature ordered phases 
are possible, and they are separated from the high-temperature disordered phase by a 
line (£) of continuous second-order phase transitions. One of them is homogeneous, here 
the non-zero averaged value of the order parameter is the same throughout the system. 
In another one, there is a periodic modulation of the order parameter along certain 
anisotropy axes. The line of transitions between these two ordered phases terminates at 
the line (£), and the point where they meet is just the LP. 

Finally, let us note that an enormous amount of papers dealing with different versions 
of " Lifshitz-like" theories appeared recently in quantum field theory (QFT), particle 
physics, gravitation, and cosmology (see [21,25-28,16,29] and papers that quote these 
references). Here most of the work is concentrated on models in rf-dimcnsional space- 
time manifolds where higher space derivatives are allowed, but higher time derivatives 
are forbidden and not generated by renormalization. In the "Euclidean picture" of (1.1), 
this corresponds to the situation when one of the "perpendicular" components of x = 
(x±,x\\) is associated with the time, while the remaining d — m — 1 components of x± 
as well as all m "parallel" components are considered as spatial coordinates. Obviously, 
the presence of anisotropy between time and space (as well as between space coordinates 
alone) breaks the Lorentz symmetry in such models. Hence the name "Lorentz violating 
(LV) theories" frequently used in the literature [21,30,7]. 

Up to few exceptions [21,25], usually one considers, as in [28,16], the case when all 
space directions scale in the same manner. In our notation this corresponds to fixing 
d — m = 1 in (1.1) and working on the line d = m + 1 of the (m, d) plane. On the other 
hand, of special interest in QFT are theories in the four-dimensional space-time, which 
corresponds to the line d = 4. The intersection point of these two lines is (3, 4) where we 
have the usual "d = 3 + 1" not forgetting about the fourth-order derivatives along the 
three spatial directions. 

Now, if we still are interested in integer values of m by keeping d — 4 but allowing 
that some spatial directions scale like time, we land at d = 2 + 2 and d = 1 + 3. Such 
possibilities have been considered in the classification of [21]. Moreover, in the special 
case m— 1, d — 4 Anselmi [25] presents explicit calculations and results for "Lifshitz type 




where A$ = - (d — m + 9m — 2 + 77^2) (1-2) 
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models" in the large-n expansion, which in the case of scalar fields completely parallel 
that of [20] and of the present paper — in the context of the condensed-matter physics 
(for an explicit comparison of calculational techniques and results see Sec. 8). 

The model (1.1) has several key features that attracted a growing interest and led 
to numerous sophisticated generalizations in different fields of QFT. First of all, the 
inclusion of fourth-order derivatives along to spatial directions raises the upper critical 
dimension of (f> 4 models d* from the usual d* = 4 at m = to d*(m) = 4 + to/2, which 
can go up to d* = 8 at to = d. This comes along with the fact that the renormalized 
values of relevant Feynman diagrams are finite up to these relatively high dimensions. 
The same can be seen as an improvement of ultraviolet convergence of Feynman integrals 
(saying nothing about technical complications!) in fixed dimensions with respect to their 
counterparts in conventional CP-like field theories. For example, the same diagrams that 
usually diverge when d > 4 are finite up to d = 5.5 when to = 3. A similar thing happens 
also to the lower critical dimension dg(m) of models with continuous 0{n) symmetry. As 
we discussed above, now we have 2 < dt(m) = 2 + m/2 < 4. 

Further, the model (1.1) can be extended in two obvious ways. It is possible to add 
terms containing two fields and derivatives higher than four in the parallel subspace. 3 
Demanding that gradient terms of highest order, say 2z, 4 scale in the same way as 
(Vj_</>) 2 yields the modified lines of critical dimensions, d*(m) = 4 + m(z — l)/z and 
d>e,z(m) = 2 + m(z — l)/z. Moreover, higher powers of fields without gradients can be 
added, and models with g s (j) s interaction terms can be considered apart from the usual 
4 one. The engineering /i dimension of the effective expansion parameter g s is [g s ]n = 
(d* z s {m) - d) (s - 2)/2 where d* z s {m) = 2s/(s - 2) + m(z - l)/z. 

Clearly, such generalizations allow to consider a very large amount of theories with 
different d, to, z, and s, which will have very different properties depending on chosen 
set of these parameters. In particular, we see that the "weighted" power counting raises 
the engineering dimensions of coupling constants at <j> s when m > and z > 1. Hence 
it is possible that certain non-renormalizable Lorentz-symmctric models can become 
rcnormalizable at the same d when the Lorentz symmetry is broken. A classification of 
LV theories involving scalar and fermion fields has been undertaken in [21] and generalized 
to gauge fields in [35,36]. 

To give some examples we mention a scalar cf) 10 model with z = 2 and m = 3 whose 
upper critical dimension d* is d* 10 (3) = 4, see [21, Eq. (2.14)], [37, Eqs. (2.4)-(2.5)]. 
Another example with d* = 4 is a <j> 6 model with z = 2 and m = 2 [21, Eq. (2.15)]. 
Reference [37] gives explicit one- loop calculations for a similar </> 6 model with z = 2, 
to = 4 and d* — 5. The last two instances belong to the class of theories describing 
the m-axial tricritical LPs in strongly anisotropic systems with short-range interactions. 
These have been considered some time ago in [38-40]. Further generalization involving 
extra long-ranged uniaxial dipolc-dipole interactions can be found in [41,42]. 

Similarly as we noted it above by discussing the LP at z = 2, the field 4> becomes 
dimcnsionless at d = de iZ (m). In this case [g s ]fj, = 2 for any to and z and independently of 



3 Early work in this direction [31-34] in condensed-matter physics did not receive much attention. 

4 Different values of z lead to different classical values of the anisotropy exponent 6q = 1/z. However, as 
we already discussed above, the critical exponent 9 differs from its classical value 9o due to the non-trivial 
rcnormalization of the interacting theory. Hence, we find it misleading to identify the derivative's power 
in the Hamiltonian with the critical exponent as it is sometimes done in the literature. 
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s. Along with the constraint d — m+1, the condition d = di :Z (m) is satisfied at z = m. In 
the four-dimensional space-time d — 3 + 1 this reduces to z = 3. Such kind of observations 
led Horava [28,43,44] (see also [16,30,45]) to his seminal formulation of the power- counting 
renormalizable quantum field theory of "gravity at a Lifshitz point" in 3 + 1 dimensions. 
However, "beyond power-counting renormalizability... explicit calculations have so far 
been extremely limited" [45]. A few examples that give some taste are [37,46,47]. 

The aim of the present paper is an explicit large-n calculation of critical exponents vl2 
and vl/± that control the behavior of correlation lengths £j_ and £m at rn-axial Lifshitz 
points. These correlation lengths are related to directions perpendicular and parallel to 
modulation axes of the periodically ordered low-temperature phase. We have £j_ ~ \t\~ UL2 
and £ii ~ l T l ^ L4 for small reduced temperature deviations r = (T — Tlp) /Tlp from the 
LP. The other exponents jl, (Xl, and Pl derived below give the thermal behavior of 
the susceptibility, specific heat, and order parameter. The present work extends that of 
[20] where the correlation critical exponents i]L2 and 77^4 have been calculated to first 
non-trivial order of the large-n expansion. 

Hopefully, the present rather technical communication along with [15,17,8,20,6,48] 
could be useful for the bright community interested nowadays in Lifshitz-type theories. 

Before going to detail of our new calculations we briefly summarize some main results 
of [20]. 

2. Lifshitz point's correlation exponents at large n 

Below, as well as in determining the exponents r]L2 and 77^4 in [20], we employ the 
method of [49] . It consists in deriving the 1 jn expansions from self-consistent equations 
for full correlation functions of the interacting theory. The actual calculations are done 
in the critical (massless) theory using the dimensional regularization, see [50, Ch. 4.35, 
4.38]. Directly at the LP, the corresponding theory is given by the LGW Hamiltonian 
(cf. (1.1)) 

n L P^} = \ j d d - m r I d m z [(V r 0) 2 + (A z <p) 2 ] +^jd d x (|0| 2 ) 2 . (2.1) 

The values Tq p and Pq p , which locate the position of LP [11,51,20] and differ from zero 
beyond the dimensional regularization, are ignored in the following calculation, similarly 
as it was done in [20]. It is assumed that A cx 1/n for n — >• 00 as in the usual <fi 4 model [50, 
Ch. 2.19]. In the construction of the large-n expansion, this is not affected by the present 
modification of gradient terms. For brevity, we use the notation r = xj_, z = xu and hide 
the parameter ao by appropriate rescaling of the model (cf. [16]). Implications of using 
the massless theory are discussed in Sec. 6. The applicability of the large-n expansion to 
physically accessible LPs is discussed in [20]. 

Deviations from the LP can be effected in two ways. First of all, the temperature 
variation will lead to the perturbation of the Hamiltonian V.lp[4>} by the term 

AHrl<t>} = lJd d x\4>(x)f. (2.2) 

Its inclusion is needed in order to derive the thermal exponents vl2, ^L4, 7l, and so on, 
which is the main objective of the present paper. On the other hand, in order to describe 
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the crossover from the LP, the term 



AHM = ^jd d x(\7 z< j>f (2.3) 

has to be included into the LGW Hamiltonian. Taking it into account would open access 
to the crossover exponent tp and the modulation wave- vector exponent j3 q . This is beyond 
of the scope of the present paper. 

In the absence of perturbations, the long-wave asymptotic of the full pair correlation 
function in the momentum representation g[? s ^ (p, q) obeys the anisotropic scaling via 

G { ; s \p, q )=p- 2 +^G^\l, qp -°) = q -^G^ (2.4) 

Here only two of the three exponents are independent. The anisotropy index 6 is defined 
as 9 = (2 — r/L2)/(4 — Vla) = vla/vl2- In the large-n limit, the non-trivial correlation 
function Grl (p, q) reduces to the free-theory propagator (p, q): 

lim G { ; s \p, q ) = (p 2 + q i )- 1 . (2.5) 

Its explicit form can be read off from the first term of (2.1) after a Fourier transformation. 
The main result of [20] have been the large-n expansions 

„(!) „(!) i Ml) 

n n 2 n 

obtained for generic number of anisotropy axes < m < d in (i-dimensional space. Their 
1/n coefficients have been derived by considering the self-consistent equation for the 
function G^ (p, q ) . In a graphical form, it may be represented as 



o 



where solid lines denote the full propagators G^ (p, q). This is a symbolic picture of the 
self-consistent equation written up explicitly in (23)- (24) of [20]. Solving this equation 
for q = and for p = yielded the 1 /n expansion coefficients 

(m) 

m K d _ m f 4-(rf-m)(l + g 4 ) 1 
^- 2 d-mJ (l + g*)3 I(l, q y (2 ' 8j 

Q 

(d—m) 

41 = , / * with (2.9) 

,i4 m(m + 2) 7 (p 2 + l) 5 I(p,l) V ; 

p 

V 2 {p 2 ) = 3(8 - m)(6 - m) + 5(m 2 + 2m - 96)p 2 + (m 2 + 50m + 144)p 4 - m(m + 2)p 6 , 
and 

^)=-ryi 1 2 )/4 + 4 1 4 ) /8. (2.10) 
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The function I(p, g), 



(d—m) (m) 



p> q> 

represents the "elementary bubble" [52] at LP. It associates with the Feynman diagram 
XHIX °f the <f) 4 theory. For any D-dimensional momentum integrals we use the notation 

(D) D 

I m s / (2 - 12) 

k 

When the integrand is rotationally invariant, the angular integration is trivial and we 
have 

J /(fc 2 ) = K D f dkk D -'f(k% where K D ee ^ = 2 (2.13) 

k 

and Sx> = 2 7 r' D / 2 /r(L>/2) is the surface area of a D-dimensional sphere of unit radius. 

In [53] we have shown that the 1/n coefficients (2.8) and (2.9) agree for arbitrary < 
m < d with the 0(e 2 ) terms of the epsilon expansion of i]L2 and i]L4 obtained previously 
in [15] and [17]. Let us recall that in the LP theory the upper critical dimension is given 
by the line d*(m) = 4 + m/2 and e = d*(m) — d. 



3. Identification of the critical exponent vl 2 



Up to now, we discussed the theory (2.1) exactly at the LP. Further, we consider a 
small temperature deviation from criticality. Following [49], we take into account the 
term (2.2) by treating it as a small perturbation. The classical momentum dimension of 
the variable t is [r] = 2. Its full scaling dimension 5 is 1/vl2- In the spherical limit [9], 
= d — m/2 — 2. Our aim will be to calculate the 1/n correction to this value. 

Let us consider the theory involving the "critical" LGW Hamiltonian (2.1) and its 
perturbation (2.2). In the presence of AH T , the full propagator gets the new tem- 
perature variable r. Now, it obeys anisotropic scaling via 

G { ; s \p,q;T) =b 2K *Gf\bp,b e q -Tb 1 /—) with 2A = 2-7? i2 . (3.1) 

~ fas) 

This implies that the generalized homogeneous function G^ can be expressed in differ- 
ent scaling representations. In particular, for b = p^ 1 , 

G { ; s \p, q ;r)=p-^G { ; s \l, q p- e ;rp-^). (3.2) 

Alternatively, the powers q^ A+r > L4 (cf. (2.4)) or t~ 71 - can be scaled out, accompanied 
by appropriate scaling functions. Following the approach of [49], we have to consider 

5 For a complete renormalization-group analysis of the critical behavior at LP see [15]. 
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infinitesimal temperature deviations form the critical theory and follow the changes they 
produce in the self-consistent equations. 

As usual, it is expected that the asymptotic small-r expansion of (p, q; r) at fixed 
finite momenta p and q contains both integer and non-integer powers of r. Including only 
the first few terms, this expansion can be written as 

G { ; s) (p,q;r)^p- 2 ^x (3.3) 
D (qp- e ) + Tp-V^'D^qp-*) + T l - aL p-^- aL) ^ L 'D(qp- e ) + 0(t 2 ) . 

The specific heat exponent oll is [9] 



d- to/2 -4 
d - to/2 - 2 



0(0 



+ 0(n- 1 ) < 0, 



(3.4) 



and the power 1 — aj, = 2/(2 — e) + 0(n 1 ) of r in the second correction term is greater 
than 1. The powers of p are given by 

— = d-^-2 + 0(n- 1 )=2-e + 0(n- 1 ) and i^ = 2 + 0(n- 1 ), (3.5) 
vli 2 v L2 

while in the epsilon expansion both of them start with 2 + 0(e). The scaling functions 
D , Di, and D possess finite limits at q = 0, and in this case 



Gr(p,0;T)~p- 2A * A+BTp- 1 l VL2 +CT 1 - aL p- { - 1 - aL) l UL2 +0(r' 2 ) 



(3.6) 



This is in full analogy with the well-known formulas from the CP theory [54-56,13,50]. 
There, the constant coefficients A, B, and C have been calculated both in the epsilon 
[54,55] and large-n [56] expansions. We shall not quote their explicit expressions, but 
only note that for n large, B = 0(l/n) while 6 A and C = const + 0(l/n). This is a 
quite general property concerning only the n dependencies in </> 4 theories. It certainly 
holds also in the LP theory. This implies that in the large-n limit (3.6) reduces to 



lim G( as) (p>0;T)~p- 2 A co +C co r 1 -«V 2 + O(r 2 ( 1 -^); 



On the other hand, (p, q;r) at n — > oo gives the Gaussian propagator 



G%(p,q;r 00 ) = (p 2 +q 4 + r oo y 



(3.7) 



(3.8) 



where ~ r 7 ~ is the inverse susceptibility of the spherical model with [9] 7^° = 
2/ (d— rn/2 — 2) = 2/(2 — e). We see that in the large-n limit the susceptibility exponent 
7^° coincides with the power 1 — cc£ as it follows from (3.4). Therefore, the small- 
temperature expansion represented by (3.7) directly matches an analogous expansion of 
the function G^{p, q; r^) from (3.8). Indeed, hence we get 



G™ (p, 0; roo ) = p- 2 [1 - rooP - 2 + Q(rl 



(3.9) 



In Eqs. (4.108) of [50] the coefficients Ci and C3 are erroneously interchanged. 
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Since the correction terms in (3.7) and (3.9) are the same, by comparing these equations 
we see that Aoo = 1, and —Coo is actually the amplitude of the inverse susceptibility in 
the large-n limit: = —C 00 t~ ,l . 

A summary is that the correction term <~ T ^-otLp-(^-a L )/vL2 m ^j^g asymptotic ex- 
pansions (3.3) and (3.6) of the full propagator G^ directly matches in the large-n limit 
the linear r^p" 2 contribution of (3.9) stemming from the spherical model's correlation 
function G^(p, q). By contrast, the linear 0(t) correction of (3.6) disappears at n — > oo 
due to vanishing of its amplitude, and it has no counterpart in the spherical limit. 

From the above consideration we draw an important practical consequence: In the 
following we shall parametrize the function by the temperature variable r 1_QL 

instead of r. Thus we shall write 

G< as) (p, q; r) = G ( f (p, q; r) = p~ 2 ^G^ (1, qp~ e ; rp~ 2 <) (3.10) 

where 

r = T X - aL and 2C=^^. (3.11) 

Directly at the LP we have the obvious relation G^ (p, q; 0) = G^ (p, q; 0) = G^ (p, q). 
Assume (cf. (3.5)) that the exponent 2( has a large-n expansion 

2C = 2Co + — + 0(n- 2 ) with Co = l (3.12) 

n 

and unknown 1/n coefficient Ci • Its determination will be the goal of the following section. 
To get a direct relation between £i and the correlation length exponent v L2 we use the 
hyperscaling law [1] 

oil = 2 — (d — m)i>L2 — mi>LA = 2 — (d — m + 9m)uL2 (3.13) 
in (3.11). This yields 

_L = d -™-2 + m — -^ + 0(n- 2 ). (3.14) 
VL2 2 n n 

We shall need also relations providing the one-to-one correspondence between the full 
propagator and its spherical limit's counterpart. At n — > oo, the function D{y) that 
appeared for the first time in (3.3), is given by 

lim D(y)=°S^, (3.15) 

which follows from (3.8)-(3.9). The constants C and A from (3.6) arc defined as 

C = D(0) and A = D (0) = G*J s) (1, 0; 0). (3.16) 

In the following section we shall find the value of (i, and hence the critical exponent 
vli up to order 0(l/n). 
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4. Self-consistent equation and critical exponents 



In the presence of a perturbation, the self-consistent equations for full correlation 
functions are of the same general form as that in the theory at criticality [49,50]. That 
is, when the deviation AH T [<fi) from the LP is included, the self-consistent equation to 
be solved below, is again given by (2.7). The only essential difference is the presence of 
an additional parameter f as the new argument of the propagator . Explicitly we 
have 

(d—m) (m) „ / \ 

[G * (P ' 9;T)] = n J J WW^) (4J) 

p' q' 

with 

(d—m) (m) 

F(p,q;r)= J J G { ; s \\p' +pl\q' + q\-T)Gf\p'^-f) (4.2) 

p> q> 

where we use the representation (3.10) for the full correlation function G^ s \ 

Following the lines of [49], we expand the both sides of (4.1) to linear order in f. 
Equating the zeroth-order terms reproduces the self-consistent equation (23) of [20], 
which yields r]L2 and t\la to order 0(l/n). Matching the 0(f) contributions leads to the 
equation 

dfG^- 1 (p,q;r)\r=o = \ [EiM - ^E 2 {p,q)] (4.3) 



with 



(d—m) (m) ^ i \ 



and 



(d—m) (m) 
p> q> 

The function W appearing in the integrand of E 2 (p,q) is given by 

(d—m) (m) 

W(p,q)= J J G { ; s) (\p'+p\,\q' + q\)drG { ; s) (p', q ':T)\ f=o . (4.6) 

p> q> 

In (4.4) and (4.6), df means the partial derivative with respect to f. 

All functions appearing in above equations are generalized homogeneous functions. 
They can be written in the scaling representations similar to that of (2.4). Thus, using 
(2.4) and (3.10) we write W(p,q) as 

W(p, q) = /-™+e™-4A,-2C W(l,p-°q) . (4.7) 
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This will be needed below, along with scaling representation (2.4) for G^ (p, q) and [20] 

F(p, q) = F(p, q; 0) = pd-m+Bm-Ah^ F{l,p- e q) . (4.8) 

Similarly as in [20], we consider the equation (4.3) at zero external momentum q = 0. 
Then, taking into account (3.6), we obtain — p 2A *~ 2l >C/A 2 on the left. The constants C 
and A are given by (3.16). Using the scaling representations for all functions involved in 
E\(p, 0) and -E 2 (p, 0), we scale out the same power p 2A *~ 2 f on the right. Matching the 
amplitudes at p 2A *~ 2 C on both sides we get the equation 

= ^[£i(l,0)- 2252(1,0)] ■ (4.9) 



Here we have 



(d—m) „ (m) A 

r 4A*-(d-m) r £) (p 9 \p + 1\~ 6 q) 

Ei 1,0 = / — = / -MLliL — — 11 4.10 

P <? 

(d—m) _ (m) 

, p 4A,-2 C -( d -m) r p 6 q WM 

E2ih0) - J ip+ipA, J G * ^fttfJmw- (4 - u) 

p 1 

As discussed in the preceding section, at large n the amplitudes A and C are of the 
same, zeroth order in 1/n, and the whole equation (4.9) has to be of order 0(1). Hence, 
the integrals (4.10) and (4.11) must have simple poles in 1/n for the compatibility of the 
left- and right-hand sides of the matching condition (4.9). The same property has also the 
analogous CP's self-consistent equation [49] 7 derived in the momentum representation. 

Similarly as in the calculation of 77^2 in [20], the origin of the singularities of momentum 
integrals in _Ei(l,0) and £^2(1, 0) at n — > 00 is the behavior of their integrands at large 
momenta p. Indeed, as p — > 00, the arguments p e \p+l\~ e q of the scaling functions D and 
in (4.10) and (4.11) tend to their limiting value q, and thus the inner q integrals are 
p-independent in this limit. At the same time, the decoupled outer d — m-dimensional p 
integrations involve the powers p 2A «- 2 C-(ci-™) both in (4.10) and (4.11). They develop 
the ultraviolet pole singularities controlled by the small value (rf^l + Ci)/ n f° r large n. 
This is easily seen by taking into account the definitions of and ( from (3.1) and 
(3.12). The upshot is that at large n the matching condition (4.9) reduces to 

-Ceo = 2Rcs 1/n ^ [£?i(l,0) - 2£ 2 (1,0)] (4.12) 

where we used (see (3.9)) the large-n limit Aoo = 1 on the left. 

Of crucial importance here is the mutual cancellation of the 0(1) contributions from 
2A^ and 2£, such that their difference is 0(l/n). This is provided by the identification 
(3.11) of 2Q along with treating the full propagator G^ in the form (3.10). An attempt 
to use instead the more standard parametrization (3.2) for G^ and to identify 2£ with 



7 By contrast to the present work, the calculations of [49] have been done in the coordinate space. That 
is why, there is no direct one-to-one correspondence in intermediate details of both calculations, even in 
the isotropic limit m = 0. 
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l/^i2 (see (3.5)) would prevent the appearance of 1/n poles in the integrals E\(1,0) 
and ^(1,0) due to the lack of this cancellation. This is different from the situation in 
the CP theory. There [49], [50, Ch. 4.38], the equation for the "correction exponent" 
2A is symmetric under the changes 2A d — 2A. This implies that the scale dimension 
2A could be identified with both powers of correction contributions, 1/v or (1 — a) /v. 
Indeed, owing to the hyperscaling law a = 2 — dv, one has (1 — a)/u = d — 1/v. 

As it was readily mentioned, the residua of the integrals E\(l, 0) and E 2 (l, 0) at n — > 00 
depend on the sum of 1 /n coefficients rj^\ and £i . We have 



(1,0) 



nK d _ m f D(q) 



+ 0(1), (4.13) 



(m) 

_ nK d 



O(l). (4.14) 



The value ryll nas been calculated in [20] and quoted in (2.8). The new coefficient Ci can 
be determined via (4.12) using the pole parts of i?i(l,0) and £"2(1,0) from (4.13) and 
(4.14). Reducing there the scaling functions in the integrands to their spherical model's 
counterparts via (3.8) and (3.15) we obtain 

(m) (m) 

Ci - -® - 2AV ™ / m + iK - /ww (4 ' 15) 

q 1 

The function 7(1, q) = limbec F(p, q), has already been encountered in (2.11). Similarly, 
J(l,q) is defined by 

(d—m) (m) 

J(p,q) = C^ 1 lim W(p,q) = f [ . ,„ 1 ,,. „ — j^— ^. (4.16) 

00 n->oo 7 J [p n +q' A ) 2 \p' +p\ 2 + \q' + q\ A 

p' q' 

The non- universal amplitude , related to the normalization of the temperature vari- 
able f, cancels in the equation for (1 as it should. 

Let us introduce the short-hand notations for frequently appearing integrals: 

(rn) 

7 1 1 

M k = K d _ m / - — — — , k integer, >2; (4.17) 
J (1 + g 4 ) fe 1(1, q) 
q 

In terms of these integrals we can write 

Ci = -Vl2 - 2A/" 2 + 4A/"j = - + 4A/j (4.19) 
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where the second equality follows by eliminating the coefficient rf£l Yl& (2-8)- Using 
(3.14), we end up with 



m 1 r 



+ m ^ + 2A/2 - 4A/j +0(n- 2 ). (4.20) 



The 1/n coefficients 77^2 and #W are known from [20] and quoted in (2.8) and (2.10). 

This is the central result of the present paper. The knowledge of along with ijL2 and 

i]L4, allows us to derive critical exponents of the susceptibility and parallel correlation 
length, -fL and i>l4 through the scaling relations [1] 

7l = v L2 {2 -i] L2 ) = v Li {A-ri Li ). (4.21) 

From the first equality we get 

(1 - in 12 - 2 V ' d - 111 '2 - 2 ' 7 n, 



= , o ( 1 - , m /o o C 7 - ) + 0(n- 2 ), (4.22) 



where 

C^^vft + ^ft+M-Wj (4-23) 

is the 1/n coefficient of the inverse value 7z ; 1 = (d-m/2-2)/2 + C 1 /n + 0( n - 2 ). The 
critical exponent ^4 is obtained from the second equality of (4.21): 

with 

C Vi = -\{d-m-2)til+N 2 - 2Mj . (4.25) 

Analogous results for and fi L follow from (3.13) and [1] /3l = (2 — 7l — a L )/2. 

Below, our findings will be checked in two limits m — > and m — > d, and in e expansion 
to order 0(e 2 ). In the special case m = 1 we shall give analytical calculations at d = 4 
and numerical results at d — 3. 



5. Large-n expansions in isotropic limits 

In this section we consider the limits m — > and m — > <i. In both cases the system 
loses its spatial anisotropy. That's why we call them isotropic. 

5.1. Critical-point limit 



When m — > 0, our system reduces to that with isotropic short-range interactions in 
the vicinity of the usual CP. In the same manner as in [20], we obtain from (4.20) 



lim v L \ = d — 2 + n + 2Kd lim 



m->0 



1 _ 2 /(i.o) 



1(1,0) I 2 (1,0). 



+ 0(n- 2 ). 



(5.1) 
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Here r\ — rf 1 ^ /n + 0(n 2 ) is the Fisher correlation exponent at the CP. In writing it 

used the result lim m _).o rf^l = °f [20]. 

The limit m — > reduces 1(1, 0) and J(l, 0) to standard integrals of the form 



wc 



(D) 



-2a 



-26 



Vb(a,6). 



(5.2) 



In the dimensional regularization, 

V In M r (-°/ 2 - fl ) r(£>/2-b) T{a + b-D/2) 

Thus, 7(1, 0) and J(l, 0) arc identified as Vd(l, 1) and Vd(2, 1), and we get 



lim i/ L 2 = = d - 2 + 2 

m— >0 



( rf -l)( rf -2) (1) l +Q(n - a)- 
4 - d n 



(5.3) 



(5.4) 



This coincides with the well-known large-n expansion result for v, the usual CP correla- 
tion length exponent: see e.g. [52] or [50, (4.294)-(4.295)]. 

Along the same lines, the CP susceptibility exponent 7 is reproduced, 



lim 7£ = 7= —— 
m->0 d—2 



1 - 



d-2 

in full agreement with [52]. 



1 



d-2 
3 d 



I 11 



V d (l,l) 



2K dT?I 



V d 2 (i,i) 



-\+0(n- 2 ) 
n 



24-d 



n 



+ 0(n- 2 ), 



(5.5) 



5.2. Isotropic LP limit 

Another non-trivial check is provided by the isotropic LP (ILP) limit, m — > d. This 
would be useless for 77^2 and vl2 separately, since each of them loses its physical signifi- 
cance at ILP. But the exponent 7l, while being combined from 77^2 and ^2 via (4.21), 
is well defined and meaningful at m — d. Its m — > d limit is expected to reproduce the 
known result of [9] (see also [57]) derived in the ILP theory, 



7ilp 



1 Chls 



d-4 n(d/4-l) 2 
The coefficient 



(10-d)(d-5) 3(d-8)(d-6) 



4(d+2) 



+ 0(n- 2 ). (5.6) 



= r(d-4) 
HLS ~ r(d/2)r 2 (d/2 - 2)r(4 - d/2) y ' ' 

coincides with that in front of square brackets in Eq. (3) of [9]. Besides the simple explicit 
factors in (4.22), we must produce the non-trivial limit m — > d of the combination C 7 
from (4.23). This has to be compared with the coefficient at — (d/4— l)~ 2 /n in (5.6). 
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The first term in C 7 disappears due to the factor d — m, while 77^2 remains finite in 
the limit m — » d 8 . It is straightforward to deal also with the second term of (4.23), since 
we have already worked out the limit m — > d for 77^4 in [20]. We have 

,. 4 1 ] 3 (rf-8)(d-6) 

^ d — = 4 —2 ° HLS - (5 ' 8) 

The explicit d-dependent factor here matches the third term in square brackets of (5.6). 

The third term in (4.23) is the integral N 2 from (4.17). The function 7(1, q) in its 
integrand has the scaling property [20] 

I(l,q) = q- 2£ I(q- 2 ,l) (5.9) 

for generic m and d, while e = 4 + to/2 — d (note that e is not supposed to be small in 
this section). With this in mind we write 

00 

/ -2(d-m)-l 1 



For vanishing d — m and q — > 00, the factor qi -2 ^ - " 1 ) -1 leads to a logarithmic singularity 
of the last integral, while the remaining part of its integrand remains regular. Hence 

lim AA = lim Kd- m K m — -— = — - = Crls , (5.11) 

m ^d m ^d drn m 2(d-m) 7(0,1) 2V d (2,2) 6 ' 1 ' 

where the explicit value of Vd(2, 2) is given by (5.3) 9 . Thus we reproduce the first con- 
tribution, 1, in square brackets of (5.6). 

The treatment of the integral Nj is similar. Here we have to take into account that 

J(l,q) = q- i -^J(q-\l) (5.12) 

where the function J(q~ 2 , 1) has a regular limit as q — > 00, namely V^(4, 2). Hence, 

lim -2AG = -K d ^H =\{10- d)(d - h)C HLS (5.13) 

m^d [Vd(2,2)] 3 

in agreement with the last remaining contribution of (5.6). 

Thus, the large-n expansion of the susceptibility exponent 7l of the anisotropic m- 
axial LP theory (4.22)-(4.23) possesses correct limits of the fully isotropic LP and CP 
for generic values of d in both cases. 

From (4.24)-(4.25) we also derive the ILP limit of the parallel correlation length ex- 
ponent UL4- 

1 1 ±C H ls 



^ILP 



d-4 n(d-4) 2 



1 | (10-d)(d-5) | 3 (d-8)(d-6) 



O(n"0. 



d(d + 2) 

The expression in square brackets differs from that of (5.6) only by the factor d in place 
of 4 in the last denominator. 



8 The same holds for 77^2 m the e expansion [17] . 

9 The same result can also be reached by considering the limit m — > d of the explicit expression for 
1(0, q) derived in Appendix B of [20] for generic m and d. 
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In the following section we shall verify our results without switching off the spatial 
anisotropy, by checking them against the existing e expansions for m-axial LPs. 



6. Compatibility of 1/n and e expansions 

In this section we show the equivalence of 1/n and e expansions for in their 
common region of validity. The compatibility of these two types of expansions for other 
LP exponents then follows from scaling laws (see e.g. [15]) by taking into account that 
for correlation exponents 77^2 and 77^4 this property has already been proved in [53]. 
First we give a rather formal development dealing with generic integral representations 
of involved functions valid for arbitrary m € [0, d]. In the following subsection we provide 
an example of an analytical calculation at to = 2 where all these functions are taken in 
their explicit forms. 

6.1. Generic to 



The critical exponent vj^ was obtained to 0(e 2 ) in [17] using dimensional regulariza- 
tion and minimal subtractions of e poles for n € [0, 00 [ and to e [0,d\- For large n it 
reduces to 



-In 6 £ 2 

L2 n 2n 



to — 4 , . TV (to) , , n 



e A + 0(e'\n- A ). 



The precise expressions for the functions j<p(m), j a (m), and J u (m) in terms of single 
integrals, as well as their explicit values at m — 2 arc listed in Section 3.2 of [17]. 
The same result can be expressed in terms of correlation exponents as 

vl\ = 2 - £ + 6- + f] L2 + m(9 -i) -UJ u {m)— +0(e 3 ,n- 2 ), (6.1) 

which has the structure very similar to that of (4.20). The factor (9 — 1/2) <~ e 2 /n rep- 
resents the non-classical, anomalous part of the anisotropy index 9. The term containing 
it, as well as n L2l have direct counterparts in (4.20). The contribution J u (m) stems in 
part from the non-trivial two-loop renormalization of the <p 4 coupling constant through 

the Feynman diagram ><C{^ °f the vertex function T4. Another source of its appear- 
ance is the similar Feynman graph (with glued left external lines) of the two-point vertex 
function with a <fi 2 insertion, T^ 2 ' 1 ). It is given by 

I 1 / 777 \ 

Mm) = l--C E - -V>(2 - - J +j u (m) (6.2) 

where Ce = — V'(l) — 0.5772 ... is the Euler's constant and ip(z) is the digamma function, 
the logarithmic derivative of the Gamma function (see e.g. [58,59]). The last term reads 

Ju (m) = %^ b m [ d m u $(«; m, 0)e(«; to) (6.3) 

t m,0 J 

where b m = 2 _4_m 7r _3 / 2_m / 4 as in (B.2), the value F TOj0 can be easily read off from (A.5), 
and S D = (2n) D K D = 2ir D / 2 /T(D/2) was introduced in (2.13). $(m;to,0) stands for the 
scaling function $ from (A.l) taken at e = 0, and <d(u;m) is defined in (B.1)-(B.5). 
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Comparing equations (4.20) and (6.1) for u L2 we see that their equivalence for large 
n and small e requires that the e expansion of combination 2A/" 2 — 4A/j must be given by 

2N 2 - 4AG = 6e - UJ u (m)e 2 + 0(e 3 ) . (6.4) 

In the remaining part of this section we shall show that this statement indeed holds true. 

To deal with A/2 and Afj defined in (4.17)-(4.18) we need the e expansions of the 
functions 1(1, q) and J(l,q). The integral I(p,q) from (2.11) has an ultraviolet pole, and 
its formal e expansion can be written as 

I(P, q) =p- £ ^[l + e io(q/VP)\ + 0(e) . (6.5) 

Using the e expansion of I(p, q) outlined in Appendix A we can write 

(m) 

Kd-m f 1 



■A/2 = e J (1 + g4)2 [1 - c s - ec-_\l (l,q) + 0(e 2 )] . (6.6) 

q 

The part of A/2 involving the first two constant terms from square brackets reads 



A/;« 



1- 0^-77 V- (2- 



e 2 + 0(e 3 ). (6.7) 



Employing the integral representation (A. 4) for I (l,q) and interchanging the order of 
integrations we obtain for the remaining part of A/2 

(m) _. 

^=-^fe 2 jd d x[G^( X )] 2 e-^ J ^lzl + 0(e 3). (6 . 8) 

Here we recognize that the outer q integral is the same as in the definition of the function 
H (r,z) introduced in (B.6). Hence, it can be viewed as the Fourier back-transform of 
Ho(r, z) in the subspace of r coordinates. This allows us to write 

A/f) = -^1/1 e 2 J d d x[G (0) {x)] 2 e -ilr j d d-rr VHQ{r ^ z) e ilr' +0(£ 3^ ( g Q) 

The integrals in (6.9) can be safely evaluated at the upper critical dimension d* = A+m/2. 
Here we use the scaling form (A.l) for G(x) along with (B.4) for Ho(r, z) and arrive at 

b m J d 4 -^r- 4+ ^ e~ ilr J d 4 -^r' e ilr ' J d m u $ 2 (u; m, 0)0 (u^) . (6.10) 

The structure of r and r' integrations here is just the same as that encountered in [53], 
which produces a Dirac delta function (2n) 4 ~ m ^ 2 S(r — r 1 ). Thus we get for A/f 

c_i 7 

where in writing the last equality we took into account the definition (6.3) of the function 
j„(m) taken from [17]. Together with (6.7) this gives 

A/2 = £ - J u (m) e 2 + 0(e 3 ). (6.11) 
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The treatment of the integral Mj requires additionally an s expansion of the function 
J(p,q) from (4.16) (see appendix B). At small e, this function can be written in the 
scaling form analogous to (6.5). Again, both the pole and finite parts have to be taken 
into account in the expansion 



J(p,q) 



h-i p 



-1-e 



e l + g 4 /F 



[l + ej (q/^p)] + O(e). 



5.12) 



By contrast to I(p, q), the 1/e pole of J(p, q) is a function of p and q and has an infrared 
origin. This can be seen from (4.16) by power counting. We cannot calculate neither 
*o(<?/\/p) nor Jo{q/\/p) analytically for arbitrary m. However, in the next section we 
shall write down these functions explicitly for the special case m = 2. 

According to the decomposition (B.8) of the function J(p, q) into three parts, we have 
three contributions to Mj. The first of them, that involves J^(l,q) from (B.9), is 



M x) = -#d- m c_i 



(m) 

/ (1+^ 



q*) 2 P(l,q) 



(6.13) 



The integral in M^ differs from that in Mi (sec (4.17)) only by the power of 7(1, q) This 
implies that in the e expansion, owing to the relation 



2e 



0(s% 



P{l,q) c-i 7(l,g) c 2 _ x 
Afj 1 * 1 can be simply expressed in terms of Mi = e — e 2 J u {m) + 0(e 3 ). This gives 

M^=e-2M 2 +e 2 \\^2 ^ ° E 



m 

T 



In 2 



+ 0(6*). 



(6.14) 



(6.15) 



The next contribution to Mj, Mj 2 \ includes the 0(s°) function from (B.14). Here 
only the leading e 2 /c 2 _ 1 term from l/7 2 (l,q) is needed and we have 



M?=e 2 



(m) 

K 4 - m/2 f J( 2 )(l,g) 
l + q i 



H 2 -t) 



In 2 



0(e 3 ) . 



(6.16) 



In the last term, Mf\ we use the integral representation (B.15) for J^(l,q). Changing 



(2) 

the order of integrations, as in M^> (see (6.8)), we obtain 



A/f = 



2 K&-m/2 , 



(m) 



Ae(ii;m)G (0) (i) e~ tlT 



-iqz 



+ 0(s 3 ). 



(6.17) 



The Fourier transformations of (A.l) imply that the last q integral can be written as 



(m) 



1 + q 4 



= d 



„i-m r , c(°\r,z)e ilr ' 



(6.18) 
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We use the scaling representation (A.l) for the free propagators in (6.17) and (6.18) 
and change the integration variable z via z = u^fr in the arguments of the scaling 
functions &{z/\/r) and ^(z/Vr'). Thus, up to a constant overall factor, the value of 
Afj 3 ^ is given by the following combination of integrals: 

Jd 4 -^r- 2+ ^ e- llr Jd 4 -^r' r'~ 2 e ilr 'Jd m u 0(«; m)$(u; m, 0)$(uy^; m, o) .(6.19) 

Proceeding here by analogy with [53] and (6.10)-(6.11) we get 

Mf ] =s 2 j u (m) + 0(s 3 ). (6.20) 

Let us note that the results for both combinations of integrals in (6.10) and (6.19), 
which contain different powers of dummy variables r and r' , are the consequence of the 
following general theorem, which can be proved along the lines of [53]: 

If a Z?-dimensional Fourier transformation f(x) <-> F(k) is defined via 

(D) 

f(x)= J F(k)e ikx , F(k) = J d D xf(x)e- ikx , (6.21) 



then 

(D) 



J S - a c lvs J d D xx-f 3 g(u-^e- lwx = w- a v-Pg(u^- 



(6.22) 



provided that a + (3 = D. 

The integrals over r and r' in (6.10) and (6.19) are covered by (6.22) with v = w = 1 
and ®(uy/ r/r') or ${uy/r/r') := g(u 2 r/r'), respectively. 

Returning to A/j, we collect terms given in (6.15), (6.16), and (6.20), and obtain 

Afj = -e + 3J u (m)e 2 +0(e 3 ). (6.23) 

Finally, following (6.4) we combine the values A/2 and Nj from (6.11) and (6.23) via 
2A2 — 4A/j. This really gives the combination Qe — 14J„(m)e 2 + 0(e 3 ) as required by 
(6.4) for compatibility of 1/n and e expansions for the critical exponent vl2- 
In the following section we show how it works in a closed- form calculation at m = 2. 



6.2. Explicit calculation at m = 2 



In this section we shall do analytical calculations for u~^ 2 m the case m = 2 for n — >• 00 
and e — > 0. Contact will be made with the e expansion result of [17]. Explicit expressions 
for relevant functions will be used without appealing to formal Fourier transformations 
employed in the preceding section. 

For large n, the 0(e 2 ) expression [17] for v^l at m — 2 reads 

vl\{m = 2) = 2 - e + 6^ + (^ _ 14 In ^ ^ + 0(e 3 , n' 2 ) . (6.24) 
To see agreement of this formula with the result (4.20) we shall explicitly expand it in e. 
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In the special case m — 2 and d*(m = 2) = 5, the scaling function (see (A.l)) 

reduces to an ordinary Gaussian, 



*(u;2,0) = (4tt)- 2 e-" 2 / 4 . (6.25) 
Then, it is straightforward to calculate I (l,q) from (A. 4) by using 

d 2 u$ 2 (u;2,0)(e^ u ^_ 1 ^ = I_i_^ c -|9 2 r_^ _ ( 6 26) 

Completing the remaining three-dimensional r integral in (A. 4) and combining its result 
with the constant O(e ) term from (A. 5) via (A. 6) we obtain 

i (qr) = 1 - q 2 arctan-^ - In (l + + 2 - C E + ln47r . (6.27) 

This function has to be used in the finite part of the e expansion of 7(1, q) given by (6.5) 
at m = 2. In this case c_i = l/(327r 2 ). The g- independent contribution here corrects the 
one erroneously given in Eq. (55) of [20] (note the different sign conventions in (6.5) and 
Eq. (44) of [20]). The precise value of this constant term was redundant for the analysis 
carried out in that paper and hence did not influence its results and conclusions. 
Now, using Mathematica [60] we obtain from (4.17) 

oo 

A/2(m=2) = K 3 - £ K 2 ^- J ^ [l-eio(q)+0(e 2 )} = e - e 2 In ^ + 0(e 3 ).(6.28) 

o 

Let us turn to the function J(p,q) and the integral Afj (see (4.18) and (B.7)-(B.8)). 
The first part of J(p,q), J {1) (p,q), is given by (B.9) with h-i = -1/(32tt 2 ). The function 
j( 2 \l,q) is given by (B.10) where at m = 2 the inner integral over u is Gaussian. It is 
proportional to e~ q r . The angular integration in the three-dimensional r integral gives 
(r/p) s'mpr. Numerical constants of these two integrations cancel giving 

oo 

j( 2) (p,q;m = 2) = — / dr lnr sinpr e^ +0(e). (6.29) 
P J 

o 

The remaining radial integral can be found, for example, in [61, 865.911]. We obtain 



jW( Pi q) = J^- «_ arcta n4-^ln(p 2 + g 4 )-C B +0(e). (6.30) 
p z + g 4 L p q z 2 J 

In j( 3 \p, q), we perform the both angular integrations of (B.15) explicitly and use, along 
with (6.25), the closed form of the function 9(u; m=2) known from [17], 

2 2 2 

6(m;2) = -2[C s + ln— + E^— J J = -2Ein(^— J. (6.31) 

Here £1(2;) is the exponential integral [58] and Ein(x) is the entire function defined by 
the expression in square brackets [58]. The radial r integral in (B.15) then reduces to 



q 2 , P 1 



OO 

j , . ^ 1 ___uV 



dr sin prj (qu^/r) = - cos (6.32) 

p vp 
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as can be found out from [62, 2.12.18.7]. Thus we remain with 

oo 

1 1 f a 2 x 

J (3) (p,q;m = 2) = -^2^2 J dxe ~ X cos ~ \° E + lnx + E ^ x ) 



(6.33) 



where we changed the variable u in favor of x = u 2 /4. It is convenient here to do the 
integration by parts remembering that the derivative of E\ (x) is — e~ x jx. The calculation 
can be completed with the help of Mathematica [60] yielding the result 



h-! \q 2 



p 2 + q 4 



arctan 
p 2p p 



q 2 q 2 q 2 1 p 2 - 

arctan In — — 

2 4p 2 



0(e). (6.34) 



p A 4p' + q** 

Adding up all three contributions (B.9), (6.30), and (6.34) to J(p,q) and expressing the 
result in the form (6.12) we derive the function jo(q)- 



jo(q) = q 2 (J^ — 2 arctan q 2 + arctan — In 



1 + ^ 4 l^.i 1— 

- -Ce + hi V7T. 



\/4T^ 



5.35) 



Finally, with the help of Mathematica [60], from (4.18), (6.5), and (6.12) we obtain 



Mj(m = 2) 



-K 3 _ S K 2 - 



qdq 



C-!j (1 + 9 4 ) 2 


= -e + 3e 2 ln^ + 0(e 3 ). 



-+jo(«) + 2i («) 



+ 0(e 3 ) 



(6.36) 



Inserting into (4.20) the known values [20] of r]L2 and 77^4 at m — 2 along with the 
e expansions (6.28) and (6.36) for J\f 2 and Afj we derive the correct result (6.24) for 
1^2 (m = 2) at large n and small e. 

7. A note on dimensional regularization and analytical continuation 



Let us consider the validity region of our results. These have been derived directly at 
LP within a " massless" theory employing the dimensional regularization. This is precisely 
the way of producing large-n expansions accepted in [49]. Above we have studied three 
situations where dimensional regularization works well and leads to meaningful results. 
These included calculations of critical exponents at CP, at isotropic LP, and the epsilon 
expansion about the upper critical dimension for m-axial anisotropic LP. The results in 
all these special cases have been obtained before by other means, but also mainly with 
the help of dimensional regularization. 

As usual, the dimensional regularization implies a need in certain analytic continua- 
tions (see e.g. [12, Ch. 4]). In the theory of isotropic CPs it is possible to obtain most of 
results, at least in lower-order approximations, in closed explicit forms. Very often, albeit 
not necessarily, they contain Euler Gamma functions with negative arguments, which re- 
quire analytical continuation. This is provided by the Cauchy-Saalschiitz formula (see 
[63, Ch. 12], [64, Ch. 1]) 



T(a) 



dtt° 



-E 



1=0 



(7.1) 
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where the Euler integral with subtractions converges for — fc— 1 < Re a < —k. In practice, 
the relation T(a) = T(a + l)/a for —1 < a < is (repeatedly) used, which gives meaning 
to such Gamma functions. For example, the coefficient j/ 1 ) from Eq. (28) of [49] contains 
the function a(2) = T(d/2 — 2) where d/2 — 2 < for d < 4. The last procedure gives 
here T(d/2 — 1) where the argument is already positive for physical range d > 2. 

Accepting the approach of [49] we worked directly at LP and used the dimensional 
rcgularization. This implied vanishing of coefficients at <j> 2 and (W z <j>) 2 in the Hamiltonian 
(see (1.1), (2.1)). No subtractions like that in (7.1) appeared in the resulting integrals, and 
it was possible to use the same transformation in isotropic limits and for infinitesimally 
small values of e. 

However, generic results contain integrals normally not expressible in a closed form. 
We have already had a hint that some analytic continuations will be needed beyond 
the special cases just discussed. Indeed, working in the s expansion we saw that the 
integral J(p, q) from (4.16) exhibits an infrared 1/e pole. However, while the dimensional 
regularization was employed in constructing the e expansion, there was no problem in 
combining it with the ultraviolet 1/e pole of I(p,q) and arriving at correct final results 
for d<d*. 

As usual, infrared problems become more severe when one moves towards lower space 
dimensions. For example, it is not possible to calculate numerically Afj directly at d = 4 
or 3 and m = 1 due to infrared divergence of the integral J(l,q). In order to obtain 
sensible results the definition (4.16) of the function J(p, q) has to be extended so that it 
becomes well-defined in the whole region of interest, that is in the interval < e < 2. 

Realizing that infrared problems are typical for massless theories we have done [65] 
another calculation of 71, employing the massive theory. We used the method of Ma [52] 
where in constructing the large-n expansion the calculations are carried out not directly at 
the transition point, but away from it. An asymptotic analysis of the theory is performed 
for infinitesimally small but not vanishing "mass", the role of which plays the inverse 
susceptibility. In this way we obtained the results for LP's critical exponents of the same 
general form as above. An essential modification occurs only in the integral representation 
of the function J(p,q). The massive theory generates appropriate subtractions, similar 
to that of (7.1), in the integrand of J(p,q). This provides the analytic continuation for 
this function of the type discussed in [12, Ch. 4.2]. In the region < e < 1 it becomes, 
instead of (4.16), 



(d—m) (m) 



J (ls \p,q) 



p> q> 



(p> 2 + q' 4 ) 2 



(p' + p) 2 + (q' + q) 4 p 2 + q 4 



(7.2) 



For larger values of e, the subtraction of l/(p 2 + q 4 ) is not sufficient and the integral (7.2) 
again diverges at small integration momenta. In the complementary region 1 < e < 2 an 
additional subtraction appears, and we must use the formula 



{d—m) (m) 

?' 2 S 2 (P,9) 



/ /TO 



_(p'+p) 2 + (q'+q) 4 p 2 +q 4 

P' Q' 

where the function S2 (p, q) is defined by 



(7.3) 
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s ^ q) = l^WT^f [ (6 - m)q4 - {2 + m)p2 ^ ■ (7 ' 4) 

Conversely, this doubly subtracted form cannot be used for e < 1 since there the last 
term of (7.3) will cause an ultraviolet divergence. 

Thus, for non-infinitesimal values of e from the range < e < 2, the function J(p, q) 
has to be replaced in (4.18)-(4.25) by its appropriate "subtracted" version J^ ls ^(p,q) or 
J (p> l) from (7.2) or (7.3). Hence follow the results for critical exponents vl2, ^L4 and 
7l , which are valid, as in [20] , for generic m-axial LPs with < m < d in the whole stripe 
2 + m/2 < d < 4 + to/2 between the lines of the lower and upper critical dimensions. 

Using the analytic continuation of J(l, q) defined by (7.3) and (7.4) we have calculated 
the 1/n coefficients for uniaxial systems (m = 1) directly at d = 4 and d = 3 where e = 
1/2 and e = 3/2, respectively. These results will be discussed in the following sections. 



8. Special case d — 4. rn — 1 



In this section we consider the special case of four-dimensional systems with uniaxial 
anisotropy. As discussed in the Introduction, this can be related to quantum field theories 
with violation of the Lorentz invariance. Fortunately, when d — 4 and m = 1 all calcula- 
tions for critical indices to order 0(l/n) can be done analytically. Details of calculations 
and results can be compared with that of Ansclmi [25] . 

Here we have to use (7.2) for the function J(l,q) in the integrand of AO in (4.18). The 
integral over p' is now three-dimensional and it involves 



(3) 

/ 



1 1 



(p'2 + 5/4)2 (p, + p) 2 + {g , + q) l 8n q ,2 [ql2 + {ql + q)2] 2 + p2 ■ 



For (7.2) we get 



OO s 



c 2 + (x + qYY + p 2 P 2 



The subtraction in curly brackets provides the convergence at x 
transformations reduce the last integral to 



(8.2) 

0. Simple algebraic 



J (ls) (p,q) = ~, 



Q .a , ..9 Jo(p,q) - o J2[p,q) - 



q^+p* 



q^+p- 



Jeip, q) 



where Jk(p, q) with k = 0, 2, 6 are given by 

00 t 

/x dx 
A A with A± = 



{x±qf 



(8.3) 



(8.4) 



The integrals Jk{p, q) can be calculated with the help of the residue calculus. For instance, 
3tt _ t q 4 + q 2 ^/4p 2 + q' i -2p 2 



Jo(p,q) 



16V2 (2g 4 -p 2 )(p 2 +? 4 )VV+? 



5) 
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For the whole combination (8.3) the result is j( ls) (p, q) =p 5 / 2 J (ls) (l, q) with 

J( .., (1 , „ i <« 4 - " + E - w gEEZ . (8 . 6) 

This function is recognized in the second term of [25, (C.2)] after removing the auxiliary 
mass m 2 introduced there artificially to avoid IR problems in evaluating the diagram 
called the scalar triangle. 

As an illustration to the discussion of Section 7, in Appendix C we rederive the result 
(8.6) in coordinate representation without any subtraction like that of (7.2). There, the 
required analytic continuation is performed by other means and involves the relation 
r(x) = T(l + x)/x for a negative x. We also believe that it is worth to give detailed 
calculations of different kinds for four-dimensional systems at LP in view of growing 
interest to such systems in quantum field theory. 



Using / _1 (l,g) = 47r\/2y v /4 + g 4 + (? 2 from [20] along with (8.6) wc: get 

*~l5s ™ d 5*b- (8 - 7) 

These agree with the terms (b) and (c) from [25, (8.3)]. Further, with rf^ = 5/(97r-\/3~) 
and <9W = -4/(27ttV3) derived in [20], we obtain from (4.19) and (4.20) 



103 , ,3 301 1 



for four-dimensional uniaxial systems. 

Both Ci/ n an d the 0(1/ n) contribution in from (8.8) do not match the value 
7 CT = — 83/(97rn\/3) from [25, (8.4)]. Such value results from the combination 7 CT = 
2[4A/j — 2Af 2 + Vl2]- Up to the sign of rf^l, it is similar to 2(i/n given by (4.19), but has 
no direct counterpart in our theory. 

9. Numerical results for uniaxial systems in d — 3 

The choice m = 1 and d — 3, corresponds to the experimentally accessible case of three- 
dimensional systems with uniaxial anisotropy. Here we have to use numerical means to 
evaluate the 1/n coefficients of the exponents 7£ an< ^ v Li, given by (4.20), (4.22)- 
(4.23) and (4.24)-(4.25), respectively. 

The integral A/2 (m=l, d=3) (see (4.17)) involves the function 1(1, q) known explicitly 
from [20], 

1(1, q) = ^(1 + <Z 4 r 1/2 (4 + q A )- 1/A K(k), (9.1) 
where K(k) is the complete elliptic integral of the first kind [66, Ch. 13], and 

k 2 = -(l- ^ I > < fc 2 < 1/2. (9.2) 

The numerical result is A/" 2 (m=l, d=3) ~ 0.249788. 
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In the calculation of Afj(m=l,d=3) (see (4.18)) we have to use, in addition, the im- 
proper integral J (2s \l,q) from (7.3)-(7.4), at e = 3/2. Numerically, using FORTRAN, 
we obtain A/j(m=l, d~3) ~ 0.1096. In fact, an involved calculation in the complex plane 
yields the explicit expression 

J (2s \l, q) = + q 4 )- 7/2 (4 + q 4 )' 3 / 4 \q 2 (9 - 46q 4 + 7q 8 ) \2E{k) - K{k) 

+ {l + q i ){l + 7q i )^A + q^K{k)] (9.3) 

where E(k) is the complete elliptic integral of the second kind. The both functions K(k) 
and E(k) are related to Gauss hypergeometric functions with the argument (9.2) via [66, 
Ch. 13], [67, Ch. 11.16] 

K{k) = ^F 1 [\,\-l-k 2 ), E{k) = ^ 2 F 1 [-\, l --l-k 2 ), \k\<l. (9.4) 

With the knowledge of (9.3) we checked the numerical value of Mj{m=l, d=3) using 
Mathematica [60]. 

Using the coefficients r]^(l,3) ~ 0.306, 6>«(1,3) ~ -0.0487 and n^(l,3) ~ 0.223 
from [20] along with above values of A/2 and Afj, we obtain 



^2(1,3) 


= 2- 1.274/n + 0(n 


-% 


^2(1.3) = 


i + 0.319/n + O(n- 2 ) 


7l(1,3) = 


= 4 - 3.161/n + 0(n" 


2 ), 


7Z 1 (1,3) = 


i +0.198/n +0(n" 2 ) 


^L4(l,3) 


= 1 - 0.734/n + 0(n 




^4 1 (1>3) = 


l + 0.734/n + O(n- 2 ) 


a L (l,3)= 


=-3 + 3.283/n+0(n~ 


% 




-- - 0.365/n+O(n- 2 ) 


Pl(1,3)-- 


= i - 0.0612/n+O(n 




/3Z 1 (1,3) = 


2 + 0.245/n+O(n~ 2 ). 



Comparing these 1/n expansions with corresponding e expansions [17] we see that 
already zeroth-order terms of both, except for are quite different. Indeed, the e 
expansions of vl2, 7l, ^L4, «l and /?l start with mean-field values 1/2, 1, 1/4, 0, 1/2, 
respectively. These do not depend on the number of anisotropy axes m and actually 
coincide with mean-field values of usual CP exponents. By contrast, the large-n limits of 
these exponents are built on the basis of the combination d~ m/2 — 2 = 2 — e, twice the 
classical dimension of the field <p{x). This value essentially depends both on m and d. 
Otherwise, the values of coefficients of the 1/n corrections are quite reasonable. Except 
for a/, and a^ 1 , they are never larger compared to the corresponding zeroth-order terms. 
Hence, there is no hope that numerical estimates derived from the short series expansions 
in 1/n will be in tight agreement with that from the e expansion. Especially this refers 
to the exponents whose zeroth orders strongly differ in both these approaches. This is 
anyway not surprising. It is well-known that the truncated large-n expansions do not 
give very good numerical estimates. 

For example, what we can get from the large-n expansions of ^2(1, 3) and v~[ 2 (l,3) 
at n = 3, is ^2(1, 3; 3) ~ 1.6 as compared to the second-order e expansion's value [17] 
0.8. In the case of /3l, where the both types of series have equal leading terms, we obtain 
3; 3) ~ 0.48, while the second-order e expansion gives the value 0.3. 
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10. Concluding remarks 

The present article represents an extension of our previous work [20] on the large-n 
expansion for m-axial LPs. In addition to the 0(l/n) values of the correlation exponents 
r]L2 and r/£4 calculated in that reference, here we derived the non-trivial 1/n corrections 
of the thermal exponents vl2, ^L4, 7l, «l, and /3l- We did it by using two different 
approaches: (i) the technique of matching the asymptotic scaling forms of full propagators 
[49,50], and (ii) the method of Ma [52] of the asymptotic analysis of the correlation 
function at small but non- vanishing temperature deviations from the LP. 

In sections 3 and 4 we described in full detail the calculations employing the method 
(i). This yields the result (4.20) for the perpendicular correlation-length exponent ^L2, 
and further exponents have been obtained through the scaling laws. Starting from (4.20) 
and using the dimensional regularization it was possible to show that: (a) at to = we 
correctly recover the large-n results of the usual CP theory, (b ) the known large-n results 
for isotropic LPs are reproduced at to — d, (c) an explicit calculation at m — 2 yields the 
correct expressions of the epsilon expansion up to 0(e 2 /n), (d) the compatibility with 
the epsilon-expansion results for generic m can be established along the lines of [53]. 

However, it was still impossible to use the formula (4.20) as it stands to compute the 
values of 1/n corrections for lower dimensions d and m ^ {0, d}, which should correspond 
to different cases of anisotropic systems. Employing the method (ii) yielded subtractions 
in the integral representation of the function J(p,q), defined in (4.16), necessary for its 
analytic continuation to arbitrary dimensions d and to of the main interest: These lie in 
the (to, d) plane between the lines d = d*(m) = 4 + to/2, d = d^(m) = 2 + to/2, to = 0, 
and to = d. In this region of d and to, the relevant expressions are (7.2) and (7.3) along 
with (7.4). These, instead of (4.16), have to be used in the integrand of J\fj in (4.18) 
which, in turn, has to be inserted into (4.20). 

Proceeding in this way we performed explicit calculations in the special case d = 4, 
to = 1 considered by Anselmi [25] in the context of the Lorentz violating field theory. 
We compared our calculations and results with that of [25] in Sec. 8. 

At last, in the case of three-dimensional systems with uniaxial anisotropy (d = 3, 
to = 1), which arc of main interest for condensed-matter physics, we were able to get our 
final results for critical exponents in the numerical form. These are compared with the 
estimates stemming from the epsilon expansion. 
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Appendix A. Epsilon expansion of I(p, q) 

It is convenient to perform the e expansions of the integrals I(p, q) and J(p, q) in direct 
space. We introduce, following [15,17], the Fourier transformation of the free propagator 



26 



(2.5) via 

(d—m) (m) 

G M(x) = GM(r,z)= J j ^£= r -***( tt)l tt =-L. (A.l) 
p i 

The scaling function $(it) = $(u;m, e) is defined by the momentum integral from (A.l) 
with r — > 1 and z — > it. Its mathematical properties are discussed in detail in [17]. 
In the coordinate representation, the integral I(p,q) defined in (2.11) is given by 

I(p,q) = j d d x[G {0 \x)] 2 e-^e-ti* where J d d x = j d d - m r J d m z. (A.2) 

Adding and subtracting unity at e~ lqz and taking p = 1 we split the resulting integral 
via 7(1, q) = 7(1, 0) + 7o(l, q) into a constant part 

7(1,0) = J d d x[G (0 \x)} 2 c- llr (A.3) 

and the function 7o(l, q) vanishing at q = and finite at d = d*{m): 

I {l,q) = J d d x[G w (x)} 2 e- llr (e"^ 2 -l) . (A.4) 

The constant 7(1,0) = F m , £ /e (see [17]) has been calculated for arbitrary m before 
[15,17,20], 

1(1, 0) = (4.)-f i r(|) ^| = ^ [1 + coe + <V)] • (A.5) 

The m-dependent pole coefficient c_i = F mi o appears in Eq. (6.5). The finite-part con- 
stant Co = (2 — Ce + ln47r)/2 is the same for any m. We cannot calculate the function 
7o(l, q) in closed form for arbitrary m. Nevertheless, a formal e expansion 

7(1, q) ^ £-i{l + [ Co + c-_\h{l,q)] e + 0{e 2 )) (A.6) 

appears to be useful in our calculations of Section 6. 

Appendix B. Epsilon expansion of J(p, q) 

To calculate J(p,q), we introduce the Fourier transformation of (p 2 + q 4 )~ 2 via 

(d—m) (m) 

j {p 2 + g 4)2 = r s *(«) ■ (B.l) 

P Q 

In direct space, the function 77(x) represents a convolution of two free propagators, a 
quite complicated mathematical object. From [17] we know that the scaling function 
^(u) has the Taylor series expansion 
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with b m = 2~ 4 ~ m n~( 6+m ^ 4 . The Laurent expansion of 4' (it) is implemented by splitting 
out the k = term of the sum. At small e, it contains an 1/e pole with the residuum 
h-i = — 26 m /r(l/2 + m/4). For arbitrary m we have = — c_i where c_i is the pole 
coefficient from (A. 5). The remaining series with fc > f is of order 0(1 ). 

Hence, the function H(x) from (B.l) can be written as H(x) = H(r, 0) + H o (r, z) where 

H(r,0) = r e *(0) = h_i + C E /2 + In y^r + lnr) + 0(e), (B.3) 
flo(r, 2) = r £ [*(u) - *(0)] = 6 m 0(u; m) + 0(e). (B.4) 
The function 6(w; to) is given by 10 

_L L ^> (-—) =lim ^-'(^(u)- 2J - . (B.5) 

On the other hand, we can handle the integral from (B.l) as in Appendix A, by adding 
and subtracting 1 at e lqz . This gives 

(d—m) (m) (d—m) (m) 

/f e ipr f f c iqz —1 

J (p2 + g4)2 ^nd ff (r,*) = J c'P- y (p2 + g4)2 - (B.6) 

p <? p <? 

Let us return to the function J(p, g). In the coordinate representation it reads 

J(p,q) = J d d xH(x)G (0) (x) e-^e-ti* . (B.7) 

In order to proceed we use here the e expansion of H(x) given by (B.3)-(B.4). This yields 
three different terms contributing to J(p,q): 

J(p, q) = J« (p, q) + (p, g) + ,/(3) (p, g ). (B.8) 

The first of them, involving only the constant contributions from H(x), is trivial be- 
cause the integration over x in (B.7) gives in this case simply the inverse Fourier trans- 
formation for G(°\x). Thus, by (A.l) and (B.3) we obtain 

jW(p,q) = hJ- + C E /2 + \n^)^— i +0{e), (B.9) 
V e / p + q 

valid for arbitrary m. This is the only contribution to J(p,q) that contains a pole term. 
The remaining parts of J can be evaluated directly at the upper critical dimension. 
Taking the lnr term from (B.3) we have 

J (2) (p,<?) = h-! Jd d '- m r r - 2+m/2 lnr e~ ipr JdTu $(u) e~ iqu ^ +0{e) . (B.10) 

Here, the inner integral over it represents a Fourier transformation inverse to that given 
in Eq. (14) of [15], which is a direct consequence of the definition (A.l): We have 

y^$ (u) e-^ = _01^ (3 2 r) with M ^^_ 1 = 1 _™_£ 



(B.H) 



10 This is the correct version of the misprinted formula (D5) of [17]. 
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Hence, after angular integration in the r subspace we can write J^(p, q) as 

oo 

J< 2 >(p, ? ) = fc_ip-V J drrlnr J^pr) K^{q 2 r) + O(e) . (B.12) 
o 

In the last two equations J^x) and K^{x) are the Bessel functions of the first and second 
kind [58], respectively. Noting that the logarithm in the integrand can be represented as 
Ymv a ^(r 0! — l)/a wc can do the last integral by using [62, 2.16.21.1]. This reference 
tells us that the result for the integral J a (p, q) involving r a is proportional to the Gauss 
hypergeometric function 2-Fi(l+a/2, [i+l+a/2; fi+1; —p 2 /q 4 )- Applying to this function 
the linear transformation 7.3.1.3 from [67] and taking e = we get 



V 2/ r(2-^J 



x(p 2 + (J 4 )- 1 -^ 2 F 1 (--,l + -;2- T ; 



TO p 



P 2 + g 4 



(B.13) 



Since one of the nominator parameters of the function 2F1 is proportional to a it is trivial 
to expand this function to first order in a. By a straightforward calculation we obtain 



J( 2 ) (p, g) = lim - \j^ (p, 9 ) - J< 2 > (p, 9) 

- - In p + g - — -j 

2 8— to p z +g 4 



p 2 + Q 4 



2-F1 



1;3 



-0(e).(B.14) 



TO p^ 

T'p 2 + g 4 



The last contribution to J(p,q) implied by (B.4) and (B.8) is 
J (3) (p, q)=b m J d d x Q{u; to)G (0) (x) e~ ipr e~ iqz +Q(e) . 



(B.15) 



It can be left in the present state for the calculation with arbitrary to, and in the special 
case to = 2, d* = 5 we give an explicit expression for J^(p, q) in (6.34). 

Appendix C. J(p, q) at d — 4 and m — 1 

We start from the integral representation (B.7) for J(p,q). As in [20], we use mathe- 
matical simplifications occurring on the line d = to + 3, to which belongs the point d = 4, 
to = 1. For d = m + 3 we have 

$d= ro +3(«) = (47r)- 2+£ c^ and * d=m+3 (u) = 2& m 7r £ /V [T(-e) - T(-e,x)} (C.l) 

with x = u 2 /4. The expression for follows from its general definition (B.2). T(— e, x) 
is the complementary incomplete Gamma function defined by the integral (see [68]) 



00 

T(a,z) = J t^e-Ut 



(C.2) 
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that converges for arbitrary a when Rez > 0. The Euler Gamma function r(— e) is defined 
through the analytic continuation of the same integral at z — (see (7.1)). 

Using the scaling forms (A.l) and (B.l) for the functions G^(x) and H(x) involved 
in (B.7) and performing there angular integrations in both the three- and m-dimensional 
sub-integrations over r and u we get 



J{p,q) 



£ 
p 



I 7T 



duu 1 £ $(m)*(u) J drr 3e/2 sm(pr)J- £ (qu^/r). 
o o 



(C.3) 



This integral representation is valid for the whole line d = m + 3 with < m < 2 and e in 
the same limits. In the special case d = 4, m = 1 we have to take here and in the scaling 
functions $ and * from (C.l) e = 1/2. At e = 1/2, the Bessel function J- e (z) <~ cosz/y/z 
and the inner integral over r does not converge. Nevertheless, it can be treated in the 
spirit of dimensional regularization in order to assign to it a sensible value. Using [62, 
2.12.18.4], the r integral from (C.3) can be written down explicitly. The result is a linear 
combination of two generalized hypergeometric functions 2F3. This can be evaluated at 
e = 1/2 yielding the needed analytical continuation for the last integral over r. A simpler 
way to do the same thing is to use the formal relation 

00 00 

J drr 1 / 2 sin(pr) cos{qu^/r) := — 2— j dx cos(px 2 ) cos(qux) 



2 2 2 2 

(2p + q 2 u 2 ) cos + (2p - q 2 u 2 ) sin Q 



4p 



4p 



(C.4) 



of the needed integral to the well-defined (for p, qu > 0) integral from [59, 2.5.22.5] 
expressed in terms of simple trigonometric functions. Using this finding we can write 



J(p,q) - p- b ' 2 (2nf/ 2 \j(w) + 2wf( W )} 
where, with w = q 2 /p, 



(C.5) 



OO 

j(w) = J du$(u)*(w)^ 



cos ■ 



sin ■ 



00 

= 2- J n- 3 J dxe-' x [r(-l/2) - T(-l/2, x)] (cos wx + sinwx). 



(C.6) 



The last integral can be done using integration by parts taking into account that r'(— 1/2, x) 
—x~ z / 2 c~ x . A straightforward calculation yields, with the help of Mathematica [60], the 
result (8.6) obtained in the main text using the momentum representation. 
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